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Abstract 

We show how the linear delta expansion, as applied to the slow-roll transition in quantum 
mechanics, can be recast in the closed time-path formalism. This results in simpler, explicit 
expressions than were obtained in the Schrodinger formulation and allows for a straightforward 
generalization to higher dimensions. Motivated by the success of the method in the quantum- 
mechanical problem, where it has been shown to give more accurate results for longer than 
existing alternatives, we apply the linear delta expansion to four-dimensional field theory. 

At small times all methods agree. At later times, the first-order linear delta expansion is 
consistently higher that Hartree-Fock, but does not show any sign of a turnover. A turnover 
emerges in second-order of the method, but the value of ($ 2 (i)) at the turnover is larger that 
that given by the Hartree-Fock approximation. Based on this calculation, and our experi- 
ence in the corresponding quantum-mechanical problem, we believe that the Hartree-Fock 
approximation does indeed underestimate the value of ($ 2 (t)} at the turnover. In subsequent 
applications of the method we hope to implement the calculation in the context of an ex- 
panding universe, following the line of earlier calculations by Boyanovsky et al., who used the 
Hartree-Fock and large-JV methods. It seems clear, however, that the method will become 
unreliable as the system enters the reheating stage. 

1 Introduction 

A period of inflation in the early universe could have the desirable consequence that a general 
initial condition will evolve towards the homogeneity, isotropy and flatness which we observe. Basic 
models require the slow evolution of a scalar field from an initial unstable vacuum state to a final 
stable state. Without knowing how to perform this inherently non-perturbative calculation exactly, 
approximation attempts must first prove themselves in the simpler situation of the quantum- 
mechanical slow roll. Though this simpler problem cannot be solved analytically, the degrees of 
freedom are sufficiently few that an exact numerical solution can be found. This allows us to test 
non-perturbative methods before proceeding to a calculation for the four-dimensional scalar field. 

The quantum-mechanical slow roll was first treated by Guth and Pi £Q, who considered 
the evolution of a Gaussian wave-packet initially centred at the top of a potential hill V = 
— \mujq 2 . Following this, the Dirac time-dependent variational method was used for a poten- 
tial V — X(q 2 — a 2 ) 2 /24, first by Cooper et al. 0, who used a Gaussian wave function ansatz, and 
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Figure 1: Complex time path. 



later by Cheetham and Copeland who included the second-order Hermite polynomial in their 
ansatz. 

The work presented here is based on an alternative variational approach, the linear delta 
expansion (LDE), recently applied 0] to the quantum mechanical slow roll. The method was 
found to reproduce the exact time dependence for longer than any of the alternative methods. 

In this paper we reformulate the LDE method in terms of a path integral rather than solving 
the Schrodinger equation with some wavefunction ansatz. Since we directly calculate expectation 
values without calculating the wavefunction, we save on calculational effort. More importantly, it 
is relatively straightforward to generalize to the generating functional formalism of quantum field 
theory in four space-time dimensions. This strategy is the same as that employed by Boyanovsky 
et al. in Ref . [3| , who were able to generalize the Hartree method of Ref. [2J . Since the LDE method 
is more successful in the quantum mechanical case, we should expect it to be more accurate when 
applied to field theory. 

We first consider the slow-roll phase transition in a one-dimensional field theory (quantum 
mechanics) with potential V = — hm&q 2 . This serves as a simple introduction to the path integral 
formulation of this problem. We then turn to a potential of the form V — X(q 2 — a 2 ) 2 / 24 where we 
outline the LDE method. Finally we demonstrate the use of this method for a four-dimensional 
scalar field undergoing an instantaneous temperature quench. 

In line with previous papers on the quantum-mechanical slow roll, we characterize the dynam- 
ical process by considering the expectation value of the field operator squared q 2 (t) (now working 
in the Heisenberg picture) with respect to an initial harmonic oscillator ground state. This is 
equivalent to the zero-temperature limit for an initial thermal distribution of states with Hamilto- 
nian H — + ^muifq 2 . We formulate the problem in this way in order to facilitate our transition 
to finite-temperature four-dimensional field theory. We have 



Green functions with respect to an initial field state at time to and a final state at time to — i/3 can 
be derived from a generating functional whose time contour c passes between these two points. 
The contour must also pass through the time t at which the q 2 (t) operator is inserted. The time 
contour typically passes from to along the real time axis in the positive direction the point t or 
beyond it. It then passes back along the real time axis to to before moving in the imaginary time 
direction to to — i/3 (see Fig. QJ. 
The generating functional is 




(1) 





(3) 
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The Lagrangian L must satisfy 



L(Re{i} = t ) = ~^q{md 2 + muj 2 )q (4) 

in order to meet the initial conditions. At later times the form of the Lagrangian may change, 
modelling some external influence on the particle. 

The field boundary conditions are fixed such that q(to) — q(to — i(3) and we derive general time 
contour ordered expectation values as follows: 



(0 | T c q(h)q(t 2 ) ■ • • | 0) 1 



Z[0] 



iSj(h) iSj(t 2 ) 



■m 



(5) 
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where we take (3 to infinity in the quantum-mechanical slow roll, but in principal could choose 
any value representing some fixed initial temperature. This we do when considering the case of 
four-dimensional field theory. The method is known as the closed time path method for studying 
real time dependent Green functions. It was first conceived by Schwinger [H] and Keldish |Jj (for 
a more recent account see jS])- 

In section 2 we outline the closed time path method and apply it to the quantum-mechanical 
model. We reproduce previous results for the inverted harmonic oscillator. In section 3 we develop 
the LDE approximation method. In section 4 we apply these techniques to four-dimensional 
scalar field theory in both first and second order. For completeness we include an Appendix on 
the derivation of the propagator D(t,t'), although this material can also be found in standard 
references (jHJ, [H])- 

In the remaining pages we shall use units where h = m = 1. 



2 Inverted harmonic oscillator 

In previous articles considering the quantum mechanical slow roll the particle begins 

at a time t = when it is described by a Gaussian wave function, centred at the top of a potential 
hill. To reproduce this situation here, we consider the particle prepared at t < in the ground 
state of an harmonic oscillator potential (corresponding to a Gaussian wave function). When t = 
we suddenly change the Hamiltonian to one with a potential hill. Subsequent real time evolution 
of the particle sees it "rolling off" the top of the hill. 

First let us consider a final potential of the form V = — ^uj 2 q 2 . In terms of the Lagrangian we 
have 

L(t) = \qK{t)q 
K(t) = -d 2 ~uo 2 {t) 

uo 2 {t) = Q(-t)u$ -@(t)uj 2 f . (6) 



e«)H I 253*2 ■ m 



where we define 

OAA — J 

Re{t} < 

To solve the field theory we begin by shifting the field variable q in Eq. Q in order to complete 
the square 

q(t)^q(t)- fdt'D(t,t')j(t'). (8) 



The propagator D must satisfy K(t)D(t, t') = S c (t, if), where the contour delta function S c (t, t') is 
defined for a test function f(t) by / At' f(t')6 c (t, t') = f(t). This results in a generating functional 
of the form 

Z\j}=Z[0]exp(-iJdt'dt" \j{t')D{t',t")j(t") \. (9) 
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Performing the functional derivatives in order to obtain (q 2 (t)} we find 



(f(t)) = (o\f(t)\o) = 



Z[0] 



Sf(t) 



iD(t,t). 



(10) 



2=0 



In the zero temperature limit, the propagator is found to have the general solution (see appendix) 



where 



[d 2 +cu 2 (t)]U ± (t)=0 
and for Re{i} < the two independent solutions are 

U ± (t) = ex.p{±iu l t}. 



(11) 
(12) 

(13) 



The dynamical information of the theory is contained purely in the t/^-functions. The problem is 
essentially reduced to solving a second order differential equation. Fixing the boundary conditions 
such that 



tf ± (0+) = ^(0-) 
d t U ± (p+) = 8^(0-) 



we find the general solution to Eq. (|12J) : 

= Q(-t)e ±iUit + 0(f) ^cosh( W/ t) ± i^- sinh(w/i)J 
Putting Eqs. ifTH)!. (fTT|> and QSJ) together we find 









H 




-4) 







[cosh(2w/t) - 1] 



(14) 
(15) 



(16) 



(17) 



This is the standard harmonic oscillator result for t < 0. For t > the expectation value begins 
to grow as the particle rolls off the top of the hill. The growth becomes exponential for large t: 

<9 2 W)-^-(^l + ^exp(2^). (18) 
This is in exact agreement with Guth and Pi |Tj after carefully comparing parameters. 

3 Linear delta expansion 

We next turn to the problem of a symmetry breaking potential described by a Lagrangian of the 
form 

L(t) = Igtf( t ) ? -Mg< 

K{t) = -d 2 t -u\t) 
uo 2 {t) = ®{-t)JI - ®{t)u) 

A(t) = e(t)A. 



(19) 



We could at this stage perform a perturbative expansion in powers of A. However, we know 
that the particle is bound by the q 4 term to a region near to q = 0. If we perturb about the 



4 



Gaussian solution for (q 2 (t)), the perturbative correction must become large so as to prevent the 
exponential increase, and the philosophy of perturbation theory therefore breaks down. 

The linear delta expansion (LDE) is a practical way of improving those aspects of a perturbative 
series which lead to its divergence ^J^J. In toy models, where exact results are achievable, the 
LDE is known to produce convergent results and to do so much faster than alternatives. See, for 
instance, ^1 E] and references therein. The LDE has also been used successfully in many other 
situations, including studies of scalar theories ^3 ■ 

In practice we substitute the Lagrangian with a new <5-Lagrangian which is the same as the 
original upon setting S equal to 1 

L —t L$ = (I — 5)L + 5L. (20) 

Here, Cq is just taken to be the quadratic part of the Lagrangian, depending on some variational 
mass fj,, 

U = -\q{d 2 t - » 2 )q. (21) 

The mass [i is treated as a constant for the purpose of performing any time integrals, and /z 2 is 
taken to be equal to —uo 2 for Re{i} < so as not to interfere with the fixed initial conditions. We 
have 

L(Re{t}<0) = -±q(d?+u;l)q (22) 



and 



L s (Re{t} > 0) = - X -q{d 2 t - f)q + S 



K-^ 2 ) 2 a 4 



q -24 q 



(23) 



Any given physical quantity is calculated as a perturbative expansion up to some given order in 5. 
We then set 5 equal to 1 and choose the value of fi according to the principle of minimal sensitivity 
(PMS). For {q 2 (t)) this is 

&1 = . (24) 
d/x 

The rationale for the PMS is that, although the exact value of the quantity in question can not 
depend on /x, the expansion will have some residual y, dependence when truncated to some finite 
order. The stationary points have a special status, in that at such points this dependence is 
locally zero. At other points, where the dependence is non-zero, there is no reason to choose one 
over another. Apart from this logical justification, it has been rigorously proved in some simple 
models that the sequence of approximations 1 provided by the PMS indeed converges (exponentially 
rapidly) to the exact answer [T2l I13| . in contrast to the perturbative expansion, where /x is fixed, 
which gives rise to an alternating divergent series. Apart from these proofs of convergence, it has 
been applied successfully, in a pragmatic way, to a large variety of problems in quantum mechanics 
and quantum field theory, both in the continuum and on the lattice [TBI 116) . 

In some problems it is unfortunately the case that there is not a unique solution to the PMS 
condition, i.e. that there are several stationary points. In that event, some element of subjective 
judgement has to be exercised, such as the width of the maximum or minimum and continuity 
with known results or expected behaviour. 

In the present problem the PMS criterion provides a different constraint on [i for each final 
time that we consider. Though /x will be different for different final times, it is not considered as 
a time dependent function in the evolution up to that final time. This is the simplest and most 
natural way to implement the LDE in a time-dependent problem. 

The propagator is given as in Eq. |1 ljl: however, the mode functions are now dependent on /i 
and satisfy 

[d 2 t + 6H K 2 - 6(tV 2 ] U ± {t) = 0, (25) 

1 It is important to note that the LDE with PMS gives a sequence of approximations, rather than a conventional 
series, in which subsequent orders merely add additional terms to the series. Instead, subsequent orders also change 
the values of earlier terms. 
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with solution 



U ± (t) = e(-t)e ±wa + 0(t) ( coshfut) ± i— sinh(ut) 



(26) 



At first order in 5, the relevant Feynman diagrams (Fig. 0) can be written out to give 





Figure 2: Contributions to (<f(t)). 



(fit)) 



iD(t,t) 
6X 



+ ^ / dt'e{t')iD 2 (t',t)iD(t',t') 

^ Jc 

<5(^-/i 2 ) f dt'e(t')iD 2 (t',t). 



(27) 



Evaluating the integrals we have (for t > 0) 

I 

,2 



i 

A 



16 2 U! 2 fi 2 



1 + K 1+ M^ 1 l<n " i,2/i/i L] 
-12 ^1 - -^j fj,tsmh(2fit) - fl + ^\ [cosh(4/i<) - 1] 



[1 - cosh(2/it)] 



| /itsinh(2^t) 



(28) 



This is a remarkably simple, explicit form for (q 2 (t)) compared with the complicated implicit 
expressions given in Ref. 0j. However, we have verified that these expressions do indeed reduce 
to Eq. (2BJl. 

To proceed, we find the optimum value of Eq. (|28|l according to the PMS criterion, Eq. (|24|) . 
The result is the curve shown in Fig. [3J We have chosen A = 0.01 and wf = w 2 , = 25A/6 (recall 
that w 2 appears with a different sign in the Lagrangian). These parameters coincide with those 
chosen in [fl [21 El ED in order that we may easily compare our results. Also shown are the exact 
result, first-order perturbation theory, and the Hartree approximation of Ref. |2J. 



First-order perturbation theory is achieved upon setting \i 



in Eq. (55)1. while the 



Hartree approximation amounts to taking fi to be a time-dependent function given by fi {t 1 ) = 
uj 2 — (X/2)iD(t' , t'). This results in a cancellation between the coupling correction and the mass 
insertion, and a self-consistent set of equations 



(fit)) 



iD(t,t) = —U-(t)U+(t) 



d 2 +u 2 {t) + -iD{t,t) 



U ± (t) = 0. 



(29) 



The LDE result is seen to track the exact result for a significantly longer time than the Hartree 
result. It then overshoots, signifying that the LDE result gives a much improved description of 
the inflationary period, but does not do so well during reheating. 
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Figure 3: Slow roll in quantum mechanics: (q 2 ^)} 1 ^ 2 versus t. The first-order LDE result is compared 
against the exact result. Also shown are the Hartree-Fock results of Ref. [2j (HF), and first-order pertur- 
bation theory (PTi). 



In quantum mechanics it is possible to go to high order in the LDE by the use of recursion 
relations. The results of this exercise were given in Ref. 4,, where the calculations were carried 
out to 0(S 7 ). It turns out that the second- and third-order calculations do not exhibit clear PMS 
points, but thereafter successive orders follow the true curve more and more accurately up to 
the turnover point, but diverge beyond that point. We can hope that in field theory in (3+1) 
dimensions, the LDE will again give a good description of the initial slow-roll process. In field 
theory, however, it is not practical to go beyond second order. 



4 Scalar field theory 

Having developed our method for quantum mechanics, it remains to see how easily it can be 
implemented for the case of field theory. We consider a single real scalar field theory with time- 
dependent Lagrangian of the form 

L(t) = y d 3 :cji$(-a?+V 2 -m 2 (t))*-A*4j (30) 

m 2 (t) = 9(-i)m 2 - Q{t)m). (31) 

With appropriate choice of the parameters, this model crudely describes a sudden temperature 
quench in which the field is driven through a phase transition at time t = 0. 
Our interest is in determining the quantity 

(l> 2 (i))^ fd 3 x{&(x,t)). (32) 
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To perform the delta expansion we again define a (5-Lagrangian by 

L->L S = (1- 5)L + 6L (33) 

with 

L = |d 3 x {i$(-d 2 + V 2 + ^)<i>}. (34) 
We replace the original Lagrangian by our 5-Lagrangian for Re{t} > 0. This gives 

L{Re{t}>0) = y*d 3 x|i$(-a 2 + V 2 + Ai 2 )* 



(35) 



2 24 

Switching to momentum space, the propagator now satisfies the relation 

K p (t)D p (t,tf) = 5 c (t,t'), (36) 

where 

K p (t) = -df-ul(t) (37) 
uj 2 p {t) = e(-t)4 <p -Q(t)uj 2 f . p , (38) 

and now 

oj 2 p = p 2 + m 2 (39) 
c4 P = M 2 -/- (40) 
The propagator has the solution (see appendix) 

iD v {t u t 2 ) = 6 c (t 1 -t2)iD>(t 1 ,t 2 ) + Qc{h-H)iD<{t u t 2 ) (41) 

where 

iD>(t u t 2 ) = ] 1 [U+(h)U-(t 2 ) + e^U;{t x )U+{t 2 )] (42) 

iD<(h,t 2 ) = 1 1 [e^'PU+MU-fo) + E/-(ii)C/+(i 2 )] . (43) 
The mode functions satisfy 

[d 2 +u 2 p (t)]U±(t)=0, (44) 

with solutions 

U p (t) = 6(-i)e ±Mi;pt + 6(*) f cosh(wf.pt) ± i^- sinhfw/.pi) ] . (45) 
4.1 First order 

The same diagrams which contributed to (q 2 (t)) in the previous section contribute to ($ 2 (£)) here. 
The essential difference from the quantum-mechanical case is that the propagators now depend 
on momentum and that any loops will involve an integration over loop momenta. The Feynman 
diagrams in Fig. give 

P 



(t> 2 (t)) = / iD p (t,t) 



S{mj-fi 2 ) / dt'e(t') / iD p (t',t) (46) 
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(cf. Eq. where we have used the notation 



d 3 p 



(47) 



The momentum integrals are divergent and must be regularized. As in Ref. [5], we assume 
a scheme which leaves the contributions from stable modes (p 2 ,k 2 > p 2 ) being negligibly small. 
The dominant growth in ($ 2 (£)) is associated with the finite contribution of the unstable modes. 
In practice this means that we may perform momentum integrals in the finite range p 2 ,k 2 < fi 2 
to achieve finite results. We simply make the replacement 



1 

2^2 



p 2 dp. 



The calculations are performed in the high-temperature limit, where /3wj ;p <C 1, so that 



coth -a; i; p/3 



In this limit the first term is found to be 



1 fj_ 

P Jp < 




cosh(2w/ ;p i) — 1] 



(48) 



(49) 



(50) 



the second is 



A 

2/3 2 



df' 



sinh [ujf- p (t — t')] 



J f; P , 



p W i;p W /;P 

cosh [ujf- p (t + 1' 



1 - 



t:p 



cosh [ojf- p (t — t')] 



and the third is 



(m 2 f 



1 

4 



y 2 ) 



2/3 



J i;p"f; P 



i\p 

j 2 
J f;p 



J f;pJ 

[cosh(2w /;k i') - 1] 
[1 — cosh(2w/ ;p £)] 



(51) 



v;p. 



w/ ;p t sinh(2cj/ ;p i) 



(52) 



In the second term, the time integral has not been performed explicitly since the result is rather 
involved. 

Finally we impose the PMS constraint at each t in order to find jj, and evaluate ($ 2 (t)) 



d(^ 2 W) 
dp 



= 0. 



(53) 



For numerical calculations the units are chosen such that Ti = c = ks = nif = 1. The 
remaining parameters are then chosen in these units to be m 2 =1, T = 1/ j3 = 4^/(6/ A) (the 
initial temperature) and A = 10~ 12 . These are chosen to coincide with those in Ref. [5]. The 
initial temperature has no particular meaning, it is simply twice the critical temperature. The 
coupling must be small for this type of model of inflation due to constraints from the spectrum of 
density fluctuations. 
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Examples of A(<& 2 )/2 as functions of /z 2 for various times are shown in Fig. 0] We observe a 
single stationary point, a maximum, which moves to the left and becomes sharper as t increases. 
The motivation for the PMS criterion is that the exact answer is independent of \x. In any finite 
order of the LDE this independence can only be achieved locally. A broad maximum indicates 
that the LDE is robust, but it becomes increasingly unreliable as the peak becomes sharper. From 
Fig. 0] we estimate that the first-order LDE can not be trusted beyond about t=ll. 




Figure 4: A{$ 2 (£))/2 versus /i 2 for t = 10.5, 11, 11.5 in first-order LDE. 

The position of the maximum versus time is shown in Fig. At small times the dominant 
part of the action is the quadratic part, and the evolution is well described by perturbation theory, 
i.e. /i 2 ~ to 2 = 1. At later times, as the fluctuations of the field grow, the quartic terms become 
more important. In the context of the LDE this is taken into account by smaller values of \j? in 
the trial Lagrangian Lq of Eq. 1)34(1 . 

The results for the evolution of the field are shown in Fig. in the restricted range of t where 
the different methods begin to diverge. Though we have no exact solution to compare with, the 
results display the same qualitative behaviour as in the quantum mechanical case studied in the 
earlier sections (Fig.OJ. 

First-order perturbation theory is achieved within the LDE framework by setting /i 2 = m 2 . 
As in the quantum-mechanical case, the Hartree result can be reproduced by considering /i to be 
a time-dependent function, this time given by 



H\i>) = mj - (A/2) / iD(t',t'). 

J P 

The resulting self-consistent set of equations are 

($ 2 W) = f iD(t,t) = \j ±-U-{t)U+{t) 

Jp P Jp i;p 



5 2 +c 2 (t) + ^ / //}</./) 



P »;p 

± 



p 



U*{t) = 0. (54) 
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All methods give almost indistinguishable results up to t ~ 9. The Hartree and LDE methods 
remain close up to the classical spinodal region (where V"($) < 0, i.e. f ($ 2 (t)) > 1 ). At 
later times the LDE method gives a larger value of ($ 2 (t)) than the Hartree method. Based on 
our experience of the quantum-mechanical case, we believe that the Hartree method turns over 
prematurely and that the LDE is closer to the exact result for longer. However, to this order it 
fails to give any indication of a turnover. As mentioned in relation to Fig. 0] the LDE becomes 
unreliable beyond t ~ 11, as the PMS peak becomes narrower. 

4.2 Second order 



To second order in the LDE there are altogether six additional graphs. These are exhibited in 
Eqs. i|55H60|) . along with their analytic expressions, where the integrals along the time contour 
have not yet been performed. We have used a more compact notation for the D's, whereby £)* * 
stands for D p (t',i) and so on. 



-s 2 x 2 

4 

-S 2 X 2 
4 

-S 2 X 2 
6 



dt' / dt" I ( ^A t ' ^X 2 



c J c J P " ^/ ki J\t2 



dt' / dt" / iDiHDi'HB?' I iD(f I iDC" 

c J c J p 



dt' / dt" / l D p \D t ; t / iDtf iDi? iDp" Vl -^ 



c J c J p 



S 2 X(m 2 - fi 2 ) 



ki,k 2 



S 2 X(m 2 - p 2 ) / dt' / dt" / iD^iD^'HDp" / iD? 

J c J c J p 



(55) 
(56) 
(57) 
(58) 
(59) 



= -5 2 (m) - fi 2 ) 2 [ dt' [ dt" { iDp*iDp"*iD*'*" 

J C J C J p 



(60) 



In performing the time integrals in t", t' over the contour of Fig. 1, the result is most easily 
expressed in terms of the real and imaginary parts of the D's, or more precisely F and p, defined 
by 



F 
P 



1 



iD> -iD<) 



(61) 



In the high-temperature limit, in which we are working, the imaginary parts are much smaller 
than the real parts: 



Eitl*2 
P 



cosh(o;/ ; pti) cosh(w/ ;p t2) H ^- sinh(w/ ;p ti) sinh(o;/ ; pt2 



2 

u /;p 



(62) 



compared with 



t lt2 



w /;p 



(smh(uj f-pti) cosh(cj/ :p t2) — cosh(w/ ;p ti) sinh(w/ ;p t2)) 



(63) 
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The resulting expressions for diagrams l|55 [l -l|fi(} |l . having set 5 = 1, are: 




In most of the diagrams there are two terms involving an integral over t" up to t and another up 
to t'. In the original forms of these expressions there were severe cancellations between the two 
integrals, which made accurate integration extremely difficult. In the present, equivalent, form the 
two integrals give roughly comparable contributions, posing no difficulty for numerical integration. 

We have evaluated all the multidimensional integrals numerically, including the time integrals, 
using the NAG Fortran routine D01FCF. The most difficult diagram to evaluate is, of course, the 
"sunset" diagram of Eqs. H57[) and Q66[) . Because the integrand depends only on the magnitudes 
of the various momenta, there are two azimuthal integrations which can be trivially performed, 
leaving a seven-dimensional integral. 

The result of these calculations is that the expectation value A($ 2 )/2 now develops a PMS 
minimum as a function of ji 2 . Examples of this behaviour are given in Fig. for the same times 
as were previously shown at first order. The maximum appears to be a spurious stationary point, 
with a runaway behaviour for A(<f> 2 )/2. 

The trend of the minimum as a function of t is similar to that of the first-order maximum, 
decreasing slowly as t increases, as shown in Fig. El The resulting plot of A($ 2 (t))/2 versus t is 
shown in Figure El where in addition to the Hartree-Fock result we also show the result of the 
first-order large- TV calculation (with N = 1). We see that the second-order result now shows 
a turnover, but at a larger value of ($ 2 (i)) than that given by Hartree-Fock. This is the same 
feature that occurred in the quantum-mechanical problem, and we believe gives strong evidence 
that the Hartree-Fock method turns over too soon in (<3& 2 (t)). In this case, where there is no 
symmetry breaking, the large- N calculation differs from Hartree-Fock only in that the coefficient 
of iD(t,t) in Eq. (54) is reduced by a factor of 3. This means that this term takes longer to 
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Figure 7: A(l> 2 (t)}/2 versus for t = 10.5, 11, 11.5 in second-order LDE. 



become important and produce a turnover, so that the maximum value is considerably greater. 
The same feature occurs in the quantum mechanical problem, where the large- N approximation 
greatly overestimates the maximum value of (q 2 (t)). 




Figure 8: The PMS minimum /i 2 versus t in second-order LDE. 

We have seen that inclusion of the second-order diagrams leads to a turnover which does not 
occur in first order. It would be tempting to ascribe this turnover to the influence of the "sunset 
diagram", the first diagram to include the important effects of rescatteringfTT]. However, for 
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Figure 9: \{& 2 (t)) /2 versus t. The second-order LDE result is shown as a solid line. Also shown are the 
results of Ref. g| (HF), and large N. 

the present calculation it is not possible to single out this particular diagram from the others. 
Its distinctive role is rather to provide for dissipation and thermalization at later times (see e.g. 
|18l 1191 12()|F where unfortunately the LDE is unreliable. The Hartree-Fock method, which in 
our language corresponds to a time-dependent fi with a particular selection criterion, provides an 
example where a turnover is achieved without the inclusion of this diagram. 

5 Discussion 

The main motivation for this work was to expand upon the available machinery for tackling out- 
of-equilibrium problems in field theory. 

The linear delta expansion, applied to the quantum-mechanical equivalent of the slow-roll 
transition, has been shown to give a consistent improvement on other methods. However, the 
Schrodinger formulation of Ref. 0] can not immediately be generalized to field theory in higher 
dimensions. We have shown how to recast the problem in terms of the closed time-path formalism, 
which can be so generalized. This is an extension which has not been achieved in other treatments 
of the quantum-mechanical problem, with the exception of the Hartree method. 

As noted in 5 , the Hartree approximation cannot probe the non-linear regions of the potential. 
Moreover, the Hartree method is a one-off approximation, which is not capable of systematic 
improvement. To understand the later time behaviour and to probe the true vacuum, calculations 
must to go beyond Hartree. The LDE, a systematic expansion with a variational component, offers 
just this possibility, although for practical reasons, it would be extremely difficult to go beyond 
second order in quantum field theory. 

The main result of the paper is the formalism outlined in section 4, and F'ig. El which provides 
a demonstration of its use in the instantaneous quench approximation in four-dimensional field 
theory in flat space-time. The next obvious extension is to couple the field to the scale factor of 
an expanding Universe. 
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Appendix 

It is here demonstrated how to solve for the propagator iD(t\, t 2 ) in quantum mechanics. We shall 
need to impose constraints due to the commutation relations and the KMS boundary condition, 
but we begin by decomposing the propagator as 

iD(t u t 2 ) = e c {h - t 2 )iD > (t 1 ,t 2 ) + c {t 2 - t 1 )iD < {t u h) (70) 

where d c (t - t') = c dt"S c (t', t"). Since K{i)D{t,t') = 5 c (t,f), it is straightforward to demon- 
strate that 

K{t)D >{< \t,t') = Q. (71) 

We shall construct D > ( < ^ from homogeneous solutions to the quadratic operator K, i.e. functions 
which satisfy i4T(t)t/ ± (t) = 0. For t < 0, these have the solution U ± (t) = exp{±ioJit}. Thus, the 
most general form for D > ^ is 

iD^^M) =a > «>[/ + (t 1 )(7-(t 2 ) +b>^U- (h)U + (t 2 ). (72) 

Other possible combinations of [/ ± can be ruled out on imposing time translation invariance at 
early times. The parameters a > ^ < ^ and 6 > ^ < ^ are to be determined. To do this we begin by 
imposing the particle equal time commutation relation 

[q,p]=i. (73) 

We make the free field identification (T c q(ti)q(t2)) = iD{t\, t 2 ) and further that p = q. This leaves 

d t2 [iD > (t 1 ,t 2 )-iD < (t 1 ,t 2 )] tl=ta =i (74) 

which constrains the free parameters as follows 

ci< - a> +b> - b< = — . (75) 

LUi 

A further symmetry requirement at equal time is that 

iD>(t,i) =iD<(t,t) (76) 

which translates to 

a > +b > =a < +b < . (77) 
Finally we impose the KMS boundary condition 

iD < (t ,t) =iD > {t -i(3,t) (78) 

or 

a< = exp{w J /3}a > (79) 

b< = exp{-w 4 /3}6 > . (80) 

Eqs. J75J), (|77|) . (|79|l and l|80|) constitute 4 constraints on our 4 parameters. The set of equations 
is easily solved yielding 

a> = b< = — r i- (81) 

2u>i expjcjip} — 1 

a<=b> = ^- eX r P{ "f . (82) 
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We now have a general solution for the propagator at finite temperature. Taking the zero tem- 
perature limit we have 

iD>(t u t 2 ) = ^-U-(h)U + (t 2 ) (83) 

t D<(t 1 ,t 2 ) = ^-U + (t 1 )U-(t 2 ). (84) 

The field theory case is much the same, with mode functions satisfying 

K p (t 1 )D p {t 1 ,t 2 ) = S c {t 1 ,t 2 ) (85) 

and 

(T c 0_ p (ti)^p(i 2 )) - ViD p (t u t 2 ). (86) 
The solution for an initial state described by a temperature 1//3 is given in the main text. 
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